How to detect the pseudospin-| Berry phase in a photonic crystal with a Dirac 

spectrum 
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We propose a method to detect the geometric phase produced by the Dirac-type band structure of 
a triangular-lattice photonic crystal. The spectrum is known to have a conical singularity (= Dirac 
point) with a pair of nearly degenerate modes near that singularity described by a spin-^ degree of 
freedom (= pseudospin). The geometric Berry phase acquired upon rotation of the pseudospin is 
in general obscured by a large and unspecified dynamical phase. We use the analogy with graphene 
to show how complementary media can eliminate the dynamical phase. A transmission minimum 
results as a direct consequence of the geometric phase shift of 7r acquired by rotation of the pseu- 
dospin over 360° around a perpendicular axis. We support our analytical theory based on the Dirac 
equation by a numerical solution of the full Maxwell equations. 

PACS numbers: 03.65.Vf, 42.25.Bs, 42.25.Gy, 42.70.Qs 
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I. INTRODUCTION 

Geometric phases (also known as Berry phases) typ- 
ically appear in optics and quantum mechanics when a 
spin degree of freedom is transported along a closed orbit 
[l| ■ The geometric phase is given by the product of the 
enclosed solid angle and the spin, independently of the 
duration of the orbit (hence the adjective "geometric"). 

The spin is usually i in the quantum mechanical con- 
text, when the spin is the electron spin. In the optical 
context, the spin corresponds to the light polarization 
and may be either | or 1 depending on whether the pho- 
ton momentum is cycled or kept fixed Q. An early ex- 
perimental detection of the spin-1 geometric phase of a 
photon was the measurement of the rotating linear po- 
larization in a twisted optical fiber '3] . For electrons, the 
recently observed 0, [1[ anomalous quantization of Lan- 
dau levels in graphene is a direct manifestation of the 
geometric phase of it acquired by a pseudospin- 1 which 
rotates over 360° in a cyclotron orbit (since the pseu- 
dospin is tangential to the velocity). 

The graphene example is unusual because the spin- 
i that is rotating is not the true electron spin but an 
orbital degree of freedom with the same SU(2) symmetry, 
emerging from the motion of the electron in the periodic 
potential of the carbon atoms. Such a pseudospin is not 
tied to the fermionic statistics of the electrons and so it 
might also manifest itself in the bosonic optical context. 

The optical analogue of graphene is a photonic crystal 
with a two-dimensional (2D) triangular lattice structure. 
Haldane and Raghu 6] showed that a pair of almost de- 
generate Bloch waves (^1,^2) = ^ near a if-point of 
the Brillouin zone can be represented by a pseudospin, 
coupled to the orbital motion. The wave equation, 



uj — uod 
Vd 

d 



(la) 
(lb) 



is the 2D Dirac equation of a spin-^ particle with mass 



/X (nonzero if inversion symmetry is broken) 6]. The 
resulting dispersion relation, 

s 2 = kl + k 2 y + ^, (2) 

reduces to a double cone in the case fi = of a perfect 
lattice, with a degeneracy at the frequency ujd of the 
Dirac point. The slope cko/dk = vd is the frequency- 
independent group velocity. The upper cone (frequen- 
cies u> > uid) corresponds to the conduction band in 
graphene, and the lower cone (u> < u>d) to the valence 
band. Several analogies between the electronic and opti- 
cal transport properties near the Dirac point have been 
analysed @, Ufi 0] • What is missing is an optical way to 
directly observe the geometric phase due to the rotating 
pseudospin, analogous to the "smoking gun" found in the 
electronic cyclotron motion [4|, [5| . 

A direct analogy is problematic because there exists no 
optical cyclotron motion. One can imagine other ways 
to have a photon execute a closed orbit, but the large 
and unspecified dynamical phase is likely to obscure the 
geometric phase. Here we show how complementary me- 
dia [l(| can be used to eliminate the dynamical phase, 
resulting in a transmission minimum that is a direct con- 
sequence of the 7r phase shift acquired by the rotating 
pseudospin. We support our argument by an analytical 
solution of the Dirac equation and by a numerical solu- 
tion of the full Maxwell equations. 



II. CALCULATION OF THE GEOMETRIC 
PHASE 

The system that can isolate the geometric phase from 
the dynamical phase is illustrated in Fig. [Tj (lower two 
panels). It is the optical analogue of the p-n junction 
in graphene studied in Ref. [llj . In graphene, comple- 
mentary media are formed when the Fermi level crosses 
from the conduction band to the valence band [12|. For 
the optical analogue, we introduce a (smooth) step in the 
Dirac frequency at x = 0, so that ujd decreases from 
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FIG. 1: Lower two panels: Schematic layout (not to scale) of 
the triangular-lattice photonic crystal (with the cross-section 
of the dielectric rods shown in green) and plot of the corre- 
sponding profiles e(x) and ^(x). The two regions x < and 
x > form complementary media if the rescaled frequency e 
is an odd function of x while the mass term jj, is an even func- 
tion of x. The top panel shows a closed orbit in the photonic 
crystal, with the dashed lines indicating tunneling through 
the region of imaginary wave vector. 



for x < to oo^ for x > 0. The Dirac frequency can 
be changed for instance by varying the radius of the di- 
electric rods that form the photonic crystal. Unlike in 
the electronic case, a shift of is generally accompa- 
nied by a shift of , from v~jj to . The corresponding 
shift in the parameter e is from e~ to e + . We define the 
complementarity frequency uj c such that 
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<^ = -e 



(3) 
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As illustrated in Fig. [2j waves of frequency oj c have the 
same wave vector in absolute value in the two regions 
x < and x > 0, but of opposite orientation relative 
to the group velocity (since k and dtu/dk have the same 
sign for x > and opposite sign for x < 0). Dynamical 
phase shifts accumulated in the two regions thus cancel, 
leaving only the geometric phase from the rotation of the 
pseudospin. 

We calculate the geometric phase for the closed orbit 
shown in Fig. [1] (top panel) . Notice the negative refrac- 
tion [H|,[II[ at the interface x = where the orbit tunnels 
between the upper and lower cones of the dispersion re- 
lation (Klein tunneling). The component k y = q of the 
wave vector parallel to the interface is conserved (because 
of translational invariance in the y-direction) , while the 



FIG. 2: Left panel: Schematic view of the conical dispersion 
relations near a _R"-point, in the region x > (solid lines) and 
x < (dashed lines). The horizontal dotted line indicates 
the frequency ui c , given by Eq. J3), at which the two regions 
form complementary media. Right panel: Hexagonal first 
Brillouin zone of the triangular lattice. The K and _ff'-points 
are indicated by filled and open dots, respectively. Only the 
if-points are excited in the geometry of Fig. [T] 



component k x = k changes sign when x i— > — x. The orbit 
is reflected at the turning points x± by a mass term n{x). 
We require fj,(— x) = [i{x) and e(— x) — —e(x). Because 
^>(x) and a x ^{—x) are then both solutions of Eq. {T]) (for 
a given y-dcpcndcncc oc e lqy ), it follows that the transfer- 
matrix M(x, x 1 ) through the photonic crystal [defined by 
&(x) = M(x,x')^(x')] satisfies 



M(x,0)a x M(0,-x) 



(4) 



This is a generalized complementarity relation [Tj] (the 
original complementarity relation [1 01 ] would have the 
unit matrix in place of a x ) . 

A trajectory description is applicable if the variations 
of (i, lud, and vd with x are smooth on the scale of the 
wave length. The spatial derivatives in Eq. ([1]) may then 
be replaced by the local wave vector, — iV — > k (mea- 
sured relative to the K point). The solution is 



^ = C 



-1/2 



- E 

iq 



cos(0/2) 
e 1 * sin(0/2) 



(5) 



with k determined from £,fi,q through Eq. {2} and 
C = (/i + e) 2 + \k + iq\ 2 a normalization con- 
stant. The angles <j), define the Bloch vector B = 
(cos </> sin 0, sin c/> sin 0, cos 0) , representing the direction of 
the pseudospin on the Bloch sphere. The rotation of the 
Bloch vector along the closed orbit is indicated in Fig. [3] 
The geometric phase $ = 57/2 is one half the solid 
angle subtended at the origin by the rotating Bloch 
vector We distinguish three contributions to fi, a 
contribution f2_ from the trajectory in the lower cone 
of the dispersion relation (x < 0), a contribution + 
from the trajectory in the upper cone [x > 0), and a 
contribution fix from Klein tunneling between the two 
cones (through the interface x = 0, indicated by dashed 
lines). The Bloch vector that sweeps out il± is given 
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FIG. 3: Rotation of the Bloch vector B along the closed orbit 
of Fig. [1] with the corresponding points numbered. The full 
rotation sweeps out a solid angle of 2ir, producing a Berry 
phase of n. 



by B = (k,q,n)/s. It follows from k(—x) = —k(x), 
n(—x) = fi(x), s(—x) — —e(x) that S1+ = — fi_, so the 
two contributions from the upper and lower cones cancel. 

The contribution from Klein tunneling between the 
points ±5x has imaginary k = in. The sign of k is pos- 
itive when tunneling towards positive x (from the lower 
cone to the upper cone) and negative when tunneling to- 
wards negative x (from upper to lower cone) — to ensure 
a decaying wave oc e~ KX . 

The Bloch vector 




(6) 



rotates in the y — z plane from B + to B through the 
positive z-axis (tunneling from upper to lower cone) and 
back to B+ through the negative z-axis (tunneling from 
lower to upper cone). The value of B± of the Bloch 
vector at points ±5x follows from Eq. ([6]) with k = 0, 




e(±8x) 



B 4 



(7) 



The resulting 360° rotation of B in the y — z plane sweeps 
out a solid angle SIk — 27r, so that the total geometric 
phase acquired in the closed orbit of Fig. Q] is $ = n. 



III. DESTRUCTIVE INTERFERENCE OF 
PARTIAL WAVES 

The Berry phase of n suppresses the formation of a 
bound state at the complementarity frequency u> c - To 
show this, we demonstrate the destructive interference of 
partial waves that return to the point of origin after mul- 
tiple tunnel events. A more formal proof of the absence 
of a bound state at to c is given in App. lAl 

The scattering problem is illustrated in Fig. 0] Partial 
wave amplitudes are labelled A n for x < and B n for 




FIG. 4: Sequence of partial wave amplitudes A n and B„, pro- 
duced by tunnel events (black circles) at the interface x = 
between two complementary media. Panel a) shows a top 
view of the multiply scattered rays, panel b) shows a more 
abstract representation. Panels c) and d) illustrate the con- 
struction of the total transmission amplitude T and of the 
Dyson equation that it obeys. 



x > 0. The wave amplitudes at a tunnel event (black 
circle) are related by a unitary scattering matrix, 



= S 



A n 
B-n+1 



s = 



r a 

t r' 



(8) 



The phase shift of ir acquired in a single closed loop A n — » 
B n — > A n implies 

arg (t) + arg(i') = tt t' = -t*. (9) 

Unitarity of S then requires that the scattering matrix 
of a tunnel event is of the form 



S = 



r -t* 
t r* 



+ |i| 2 = l. 



(10) 



An initial wave amplitude A„ lal interferes with the 
sum A^ nal of partial wave amplitudes that return af- 
ter different sequences of tunnel events. Each sequence 
A n —>••■—> B n A n includes B n exactly once. 

We write A n nal = T'TA™ tia \ with T the total trans- 
mission amplitude from A" 1Itlal to B n and T the total 
transmission amplitude from B n to A^ nal . 

For T we can construct a Dyson equation (see Fig. 01 : 



T = t + r'Tr + r'Tt'Tr + r'T{t'T) 2 r + 
rr'T 

= t + 



l-t'T 
Similarly, we have 



r = t' + 



rr'T 
1-tT 



(11) 
(12) 

(13) 
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FIG. 5: Frequency dependent transmission probability for two 
values of L in the case fj,R = \il of complementary media. 
The solid curves show the numerical result from the Maxwell 
equations, while the dashed curves are calculated analytically 
from the Dirac equation. The vertical dashed line indicates 
the complementarity frequency u c . 



The two Dyson equations can be combined into a single 
equation for the variable £ = T/t = 'T'/t', 



1 - tt'Z ' 



(14) 



At this point we invoke the Berry phase relation ((9 
which together with unitarity implies tt' = — \t\ 2 = rr' 
1. The Dyson equation (fT4"|) then reduces to 



e = i/\t\ 2 . 



(15) 



Regardless of the ambiguity in the sign of £, we can con- 
clude that 



T'T = fVt = -1 => A™ 1 = —A™ 



iltial 



(16) 



The end result is therefore a phase shift of tt between 
A^ nal and A™ ltlal , without any change in the magnitude. 
The destructive interference of A^ nal and A™ ltlal , which 
prevents the formation of a bound state at frequency u> c , 
is a direct consequence of the phase shift of tt acquired 
in a single closed loop, even if the weight \t\ 2 of a single 
loop is small. 



IV. DETECTION OF THE DESTRUCTIVE 
INTERFERENCE 

To detect the destructive interference, we propose a 
measurement of the transmission probability T of res- 
onant tunneling of a plane wave through the photonic 
crystal |15| . If the confinement at x = ±L is strong, the 
transmission probability will have narrow resonances at 
the frequencies of the quasi-bound states. The destruc- 
tive interference at u> — lo c will produce a transmission 
minimum for any L. This is unlike usual Fabry-Perot res- 
onances, which would shift with L, so that there would 



FIG. 6: The same as Fig. [5] for the case fiR — —{il when the 
complementarity is broken by the mass term. 



not be a systematic minimum or maximum at any par- 
ticular frequency. 

For a well-developed conical band structure we take 
the TE polarization (magnetic field parallel to the di- 
electric rods). The parameters of the photonic crystal 
are summarized in a footnote [l^]. A 7% increment of 
the radius of the rods (at fixed lattice constant a) shifts 
lud and Vd by about 5% and 15%, respectively. The 
mass term at x = ±L is created by breaking the inver- 
sion symmetry through the addition of an extra rod in 
the unit cell (see Fig. [T|) . We have solved the full Maxwell 
equations with the finite-difference time-domain method 
[l7j | using the meep software package [l8[ . The wave vec- 
tor fe = (fc, q) is the displacement of the wave vector of 
an incident plane wave from the if-point at wave vector 
^(V^l). There are two inequivalent ii-points 
in the hexagonal first Brillouin zone, see Fig. O and we 
excite a single one by orienting the lattice relative to the 
incident plane wave as indicated in Fig. [1] [The angle 
of incidence is spread over a narrow interval 69 ~ 2.3° 
around 9 = arcsin(oK'j / /a; c ).] Results are shown in Fig. [5] 
(solid curves) for two values of L. 

As an independent test on the accuracy of the numer- 
ical calculations, we have also calculated analytically the 
transmission probability from the Diracjsquation ([T]) , us- 
ing the transfer matrix method of Ref. [jj . For simplicity 
we assumed in this analytical calculation an ideal cou- 
pling between the plane waves in free space and the Bloch 
waves in the photonic crystal. The analytical results are 
also plotted in Fig. [5] (dotted curves) and are found to 
agree well with the numerical results from the Maxwell 
equations. There are no adjustable parameters in this 
comparison. (The parameters lojj,vjj,ii were extracted 
independently from the band structure, calculated using 
the mpb software package (l9j.) 

We observe in Fig. [5] a transmission minimum at lo c 
that does not shift with variations of L. To test our 
interpretation of the origin of this minimum, we have 
broken the complementarity of the media by inverting 
the sign of the mass term at the left end of the crystal. 
(This can be done by inverting the position of the extra 
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rod in the unit cell.) For q = the inversion produces 
an extra phase shift of tt that switches the destructive 
interference to constructive interference — in agreement 
with the observed switch (see Fig. [5]) from a transmission 
minimum to a transmission maximum at u) c . 

V. CONCLUSION 

In conclusion, we have proposed a method to detect 
the pseudospin- ^ geometric phase produced by the Dirac 
spectrum in a photonic crystal. The dynamical phase 
can be eliminated by measuring the transmission through 
complementary media, so that only the tt geometric 
phase remains and a parameter-independent transmis- 
sion minimum results at the complementarity frequency. 
Our analysis is based on the Dirac equation, which is an 
approximate long-wave length description, but it is fully 
supported by an exact numerical solution of the Maxwell 
equations in a triangular lattice of dielectric rods. 

The experiment proposed and analysed here can be 
seen as the optical analogue of the detection of the geo- 
metric phase acquired during electronic cyclotron motion 
in graphene [1,0]. There is one fundamental difference: 
In a cyclotron orbit the tt phase shift is produced by 360° 
rotation of the pseudospin in the x — y plane of the lat- 
tice, while in our complementary media the rotation is 
in the perpendicular y — z plane. The difference shows 
up in the dependence of the geometric phase on a mass 
term \io z in the Dirac equation. A nonzero mass pushes 
the pseudospin out of the x — y plane, thereby reducing 
the enclosed solid angle and hence reducing the geomet- 
ric phase acquired during a cyclotron orbit [2fJ. In the 
complementary media the geometric phase remains equal 

to TT. 

In graphene, the suppression of the density of states at 
a p-n junction is analogous to the proximity effect in a 



normal-superconductor junction Observation of the 
optical counterpart presented in this paper would open 
up the possibility to study superconducting analogies in 
nonelectronic systems. 
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APPENDIX A: ABSENCE OF A BOUND STATE 
AT THE COMPLEMENTARITY FREQUENCY 

The demonstration of destructive interference of par- 
tial waves given in Sec. Illll cxplicitlv shows how the Berry 
phase of tt prevents the formation of a bound state at 
the complementarity frequency lu c . A more formal proof, 
that does not rely on the partial wave decomposition, is 
given here. 

We use again the property that if 1 $>(x)e lqv is a solution 
of Eq. (P) at u> = uj c , then also a x 1 &(—x)e lqy is a solution 
at the same frequency. We may therefore take even and 
odd superpositions of these two states to form new bound 
states \I/± that satisfy a x ^> ± (0) = ±\fr±(0). The photon 
flux density through the interface x = is 

^tf|(Q)(7 x tf ± (0) = ±^|tf ± (0)| 2 . (Al) 

This should vanish for a bound state, which is only pos- 
sible if ^±(0) = 0, meaning that the two regions x < 
and x > are decoupled. Any tunnel coupling between 
the two regions will result in ^±(0) ^ 0, preventing the 
formation of a bound state at lo c . 
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